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The adiabatic Aharonov-Bohm (AB) effect is a manifestation of the Berry phase acquired when 
some slow variables take a planar spin around a loop. While the effect has been observed in molecular 
spectroscopy, direct measurement of the topological phase shift in a scattering experiment has been 
elusive in the past. Here, we demonstrate an adiabatic AB effect by explicit simulation of the 
dynamics of unpolarized very slow neutrons that scatter on a long straight current-carrying wire. 

PACS numbers: 03.65.Vf, 03.75.Be 


I. INTRODUCTION 

Aharonov and Bohm |lj pointed out that a charged 
quantum particle may acquire an observable phase shift 
by circling around a completely shielded magnetic flux. 
This remarkable effect is purely non-classical as the elec¬ 
tromagnetic field vanishes at the location of the particle, 
which thereby experiences no Lorentz force. The origin 
of the phase shift is topological: it only depends on the 
winding number of the particle’s path around the mag¬ 
netic flux. 

The original Aharonov-Bohm (AB) setup for electric 
charge belongs to a larger class of topological phase ef¬ 
fects. This includes topological phase shifts arising for 
electrically neutral particles in certain electromagnetic 
configurations mug as well as in general quantum sys¬ 
tems that undergo adiabatic evolution fT. A paradig¬ 
matic example of the latter is a molecule that acquires 
an AB phase shift when it reshapes slowly around a con¬ 
ical intersection in nuclear configuration space |5]. This 
adiabatic AB effect is imprinted in the spectral prop¬ 
erties related to the pseudorotational molecular motion, 
and has been observed [6] in the metallic trimer Na 3 . It 
has further been predicted [7 in scattering-type chemi¬ 
cal reactions, such as in the hydrogen exchange reaction 
H + H 2 . Due to subtle cancellation effects, however, di¬ 
rect observation of the AB effect in molecular scattering 
has been elusive in the past mm- 

Matter waves in spatially varying electromagnetic 
fields is a tool to engineer a wide range of quantum ef¬ 
fects B3HE]. If the variation of these fields is sufficiently 
slow, the particle motion is governed by adiabatic gauge 
fields mm similar to those in molecules and has been 
proposed to give rise to AB phase shifts under certain 
conditions pa eng. Here, we develop an adiabatic AB 


* e-mail: erik.sjoqvist@kemi.uu.se 


effect for matter waves that scatter on a static inhomo¬ 
geneous magnetic field produced by a current-carrying 
wire. The setup uses the quantum properties of neutrons 
to induce the AB effect. Earlier experimental demon¬ 
strations of other topological phase effects, such as the 
Aharonov-C as her m and scalar AB phase shifts [20] . 
as well as related Berry phase induced spinor rotations 
EU|22], have shown that neutrons are suitable to realize 
AB effects due to their robust internal structure, the ease 
of detecting them with almost 100 % efficiency, and their 
weak coupling to the environment. Our setup is simi¬ 
lar to the one proposed in Ref. [23l for a nanoelectronic 
system. 


II. PHYSICAL SETUP 


Imagine a colimated beam of neutrons that scatter on 
a long straight wire of radius R , carrying an electrical 
current I w . The resulting magnetic field is given by Biot- 
Savart’s law 


B = 


27 tR 


f(r)e g , 



r > 1, 
r < 1, 


( 1 ) 


where we have assumed that the wire points along the z- 
axis and carries a uniform current density. Here, 0 is the 
polar angle in cylindrical coordinates with corresponding 
basis vector e#, r is the distance from the wire in units of 
R , and /io is permeability of vacuum. The magnetic field 
induces a local energy splitting of the neutron spin states, 
as described by the Zeeman Hamiltonian TL = — /x • B = 
• B with /x = —9.65 • 10 -27 J/T the neutron mag¬ 
netic moment and cr = ( a x ,a y ,a z ) the standard Pauli 
operators representing the neutron spin. The eigenvalues 
of the Hamiltonian are zb f(r) = =LVb/(r), which 

introduce a potential energy barrier (well) of height Vo 
(depth —Vo) a t the surface of the wire. 

The Hamiltonian relevant for the AB effect is given by 
the kinetic energy of the neutron in the xy plane, plus 
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the Hamiltonian describing the neutron spin, i.e., 


H --J^V 2 
Htot ~ 2 m Vr 


n. 


( 2 ) 


We write the total wave function as |^(r, s)) = 

¥>+( r > s ) |X+W) + <P-(r,s) |X-(0)), where r = (x,y) = 
r(cos#, sin#) and |x±(#)) = ^75 ( =F * It) + ^ ID) are the 
local spin eigenvectors with a z |t) = |t)> G z |i) = — ||). 
The two spin eigenstates are associated with the potential 
energies ±Vo/(r). Here, s = ht/(2mR 2 ) (dimensionless 
time) with m = 1.67 • 10 -27 kg being the neutron mass. 

To determine the spatial wave functions (/?±(r, s) we 
eliminate the spin degree of freedom by multiplying the 
Schrodinger equation with (x±\ yielding the coupled 
equations 


id s ip±(r,s) = 


-d 2 r - -d r 
r 


±K/(r)+ i 

+ (^2 +i^d e U T (r,s), 


<P±(r, s ) 


( 3 ) 


where n = 2 mR 2 Vo /h 2 . 

The adiabatic regime is characterized by negligible 
Born-Huang potential and off-diagonal coupling 

terms (^2 ± <%) </?=f(i*, 5 )- For l° w neu t ron velocities, 

the wire is impenetrable and the validity of the adia¬ 
batic approximation can be checked by comparing the 
size of the scalar potential ft/(r) and 1 /r 2 at the surface 
of the wire. This yields the adiabatic condition k 1. 
In this regime, it follows from Eq. ©> that the spatial 
wave packets ^±(r, s) are separately determined by the 
effective Schrodinger equations 


id s <p±(r, s) = 


—dr - d r 


~~2 (de + iu) 2 ± nf(r) 


if±{r,s) 


( 4 ) 


where a/r = ^ is the Mead-Berry vector potential point- 
ing in the 0 direction. 

The shift do do+ia is the key origin of the AB effect. 
It originates from the adiabatic elimination of the spin, 
which introduces the Mead-Berry vector potential mm 
A = i (x± \ V r \x±) = — (a/r)eo in the kinetic energy of 
the neutron. This defines the effective magnetic flux 


4> = j) A • dx = —2ira = —it (5) 

picked up by the phase when the neutron moves around 
the wire along a loop C. vanishes if the neutron does 
not encircle the wire. The vector potential thus corre¬ 
sponds to an AB flux line carrying half a flux unit (semi- 
fluxon) sitting at the conical intersection point at r = 0 . 
The flux-induced phase is a topological property as it 
depends neither on the detailed shape of the loop C nor 


on the dynamical parameters of the system, such as the 
electrical current I w and the speed of the neutron. A 
consequence of the topological nature of the phase shift 
is that for a non-uniform electrical current density, we 
can predict that there must be an odd number of points 
at which the magnetic field vanishes inside the wire as 
each conical intersection carries half a flux unit [24, 25 . 


III. RESULTS 

We numerically simulate an unpolarized neutron that 
scatters on the wire. We solve Eq. @ for the two adia¬ 
batic channels and calculate the density profile 

#( r > s ) = \ (l^+(r,s)| 2 + |</?-(r,s)| 2 ) (6) 

obtained by tracing over the spin. We choose I w = 10 
mA and R = 10 -5 m, which yields k = 29. The spatial 
degree of freedom of the neutron is modeled as an in¬ 
coming Gaussian wave packet of ~ 20 /im width, whose 
centre is moving in the xy plane straight on to the wire at 
relative velocity v xy . The wire is assumed to be made of 
an electrically conducting material with a positive neu¬ 
tron optical potential that can reflect slow neutrons. An 
appropriate choice would be copper with a positive op¬ 
tical potential corresponding to a critical neutron veloc¬ 
ity of about 5.7 m/s. We assume neutron velocities v xy 
well below this critical value, which means that we can 
take (p±(r,s) to vanish at r = 1. With the chosen ft, 
we expect the adiabatic treatment to be valid, which has 
indeed been confirmed numerically to a high degree of ac¬ 
curacy by comparing the wave packet dynamics arising 
from Eqs. Q and ©>• 

The summation-by-parts-simultaneous approximation 
term (SBP-SAT) method 26 is a time-stable well-proven 
high-order difference methodology suitable for solving 
wave-dominated phenomena m- In the present study a 
sixth-order accurate SBP-SAT approximation has been 
employed, solving Eq. with high fidelity. 

To study the adiabatic AB effect we first compare the 
dynamics with v xy = 0.02 m/s for a = ^ and a = 0 , 
where the latter case corresponds to ignoring the adia¬ 
batic flux. As seen in Fig. [l] (see also Supplementary 
Movies 1 and 2 ), the effective AB flux introduces a nodal 
line in the forward direction, which is not visible when 
a = 0. This is the key signature of the adiabatic AB 
effect related to the effective flux line carrying half a 
flux unit associated with destructive interference in the 
forward direction mm [28]. It would be directly de¬ 
tectable by measuring the neutron density profile in the 
xy plane behind the wire. We further perform simula¬ 
tions at v xy = 0.03 and 0.04 m/s, shown in Fig. [2] The 
nodal line persists, which demonstrates the nondispersive 
nature of the AB effect and opens up for the possibility 
to perform the test even if we allow for a distribution of 
incoming velocities within the adiabatic regime. 

A significant simplification of the experimental setting 
is that unpolarized neutrons can be used, due to the 
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FIG. 1. Scattered density of unpolarized neutrons moving with incoming velocity 0.02 m/s perpendicular to a very long 
straight current-carrying wire. Left panel includes the adiabatic AB effect. Right panel neglects the adiabatic AB effect. The 
key difference is the AB-induced nodal line in the forward direction (x > 0, y — 0) that is absent when neglecting the AB effect. 
The wire represented by the white region around the origin has radius 10 -5 m and carries 10 mA electrical current. 
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FIG. 2. Scattered density of unpolarized neutrons at incoming perpendicular incoming velocities 0.03 m/s (left panel) and 0.04 
m/s (right panel) with the same wire radius and electrical current as in Fig [I] The nodal line is clearly visible demonstrating 
the nondispersive nature of the AB effect. 


spatial bifurcation caused by the energy splitting of the 
two local spin eigenstates. We demonstrate this effect in 
terms of the local spin density 

cr ( r > s) = 1 (|<p+(r, s)| 2 (x+k |x+) 

+ k-(r, «)| 2 (X-klX-)) 

= \ (k+(r,s)| 2 - k-(r,s)| 2 ) e e . (7) 

The scalar part \ ^|y?+(r, s)| 2 — |y?_(r, s)| 2 ^ is shown 

Fig. [3] (see also Supplementary Movie 3). We note that 
only/x_(#)) contributes significantly to the AB phase 
shift; the reason being the attractive nature of the adi¬ 
abatic potential —Vof(r) associated with this state for 


r > 1. 

The low neutron velocities can be experimentally 
achieved by letting the neutron hit the wire at small an¬ 
gle. This can be realized for cold neutrons (speed a few 
100 m/s) initially moving parallel to a horizontal wire. 
The bending of the neutron beam in the gravitational 
field induces a small velocity component v xy towards the 
wire. For fixed incoming speed, the relative velocity v xy 
can be controlled either by moving the wire in the vertical 
direction or by tilting it slightly upwards in the direction 
of motion. 

The incoming wave packet used in the numerical simu¬ 
lations have been chosen to clearly demonstrate the adi¬ 
abatic AB effect on the spatial distribution of the scat¬ 
tered neutrons. However, to reproduce this wave packet 
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FIG. 3. Local spin density of AB-scattered unpolarized neutrons. The incoming perpendicular velocity of the neutrons 
is 0.02 m/s with the same wire radius and electrical current as in Fig [I] Due to the attractive (repulsive) scalar potential felt 
by the local spin down (up) state, the measurable AB effect is fully dominated by the spin down state. 


dynamics in an actual experiment would be challenging 
as it would require a very high angular resolution and 
precise centering of the neutron beam. The difficulty is 
associated with the small size of the wave packet that 
would decrease the number of neutrons that would hit 
the wire in the desired manner. A less demanding setting 
would require the transversal width of the wave packet to 
be much larger than the width of the wire as this would 
allow for higher neutron flux. The AB-induced destruc¬ 
tive interference effect in the forward direction would be 
intact for such a wave packet, as the AB effect is topolog¬ 
ical and therefore independent of shape and size of the 
neutron wave function. 


IV. CONCLUSIONS 

We have demonstrated an adiabatic AB effect for slow 
neutrons that scatter on a magnetic field produced by an 


electrical current through a very long wire. The mecha¬ 
nism of the effect is the Berry phase of the neutron spin 
restricted to the plane perpendicular to the wire. There¬ 
fore, the acquired phase shift is restricted to 7r, known 
to be the only possible non-trivial value of a planar spin. 
The 7r phase shift causes destructive intereference in the 
forward direction, providing an unambigous signature of 
the adiabatic AB effect in a scattering setup. We have 
further demonstrated the nondispersive nature of the ef¬ 
fect, which opens up for the possibility to observe the ef¬ 
fect for higher neutron velocities in the adiabatic regime. 
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